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Abstract 

Multiscale correlation functions in high Reynolds number exper- 
imental turbulence, numerical simulations and synthetic signals are 
investigated. Fusion Rules predictions as they arise from multiplica- 
tive, almost uncorrelated, random processes for the energy cascade are 
tested. Leading and sub-leading contribution, in the inertial range, 
can be explained as arising from a multiplicative random process for 
the energy transfer mechanisms. 

Two different predictions for correlations involving dissipative observ- 
able are also briefly discussed. 
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1 Introduction 



In stationary turbulent flows, a net flux of energy establishes in the inertial 
range, i.e. from forced scales, L, down to the dissipative scale r^. Energy 
is transferred through a statistically scaling-invariant process, which is char- 
acterized by a strongly non-Gaussian (intermittent) activity. Understanding 
the statistical properties of intermittency is one of the most challenging open 
problem of three dimensional fully developed turbulence. 

Intermittency in the inertial range is usually analyzed by means of the 
statistical properties of velocity differences, 5 r v(x) — v(x + r) — v(x). In 
particular, in the last twenty years overwhelming experimental and the- 
oretical works focused on structure functions: S p (r) = ((5 r v(x)) p ) . A wide 
agreement exists on the fact that structure functions show a scaling behav- 
ior in the limit of very high Reynolds numbers, i.e. in presence of a large 
separation between integral and dissipative scales, Ljr^ — > oo: 




The velocity fluctuations are anomalous in the sense that ((p) exponents 
do not follows the celebrated dimensional prediction made by Kolmogorov, 
C(p) = p/3- 111 fact, ((p) are observed to be a nonlinear function of p, which 
is the most important signature of the intermittent transfer of fluctuations 
from large to small scales. 

In order to better characterize the transfer mechanism, it is natural to look 
also at correlations among velocity fluctuations at different scales and at 
different times. The prototype of such a class of correlations is: 

C M (r, R- r) = ((5 r v{x, t)Y(6 R v{x, t + r)Y) (2) 

where hereafter we will always assume the obvious notation: r < R. 
Unfortunately, the non-trivial time dependency of correlations such as (0) is 
completely hidden, in a Eulerian reference-frame, by the sweeping of small 
scales by large scales. The "positive" side of sweeping is connected to the 
Taylor hypothesis, i.e. to the possibility of identifying single point measure- 
ments at a time delay r with single time measurements at separation scale 
r ~ rV where V is the large scale sweeping velocity. The "negative" side 
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of sweeping is connected to the fact that the inertial time-scales are always 
sub- dominant with respect to the sweeping time. This implies that in order 
to measure the temporal properties of the inertial range energy-transfer it is 
necessary to abandon the usual Eulerian reference-frame and to move in a 
Lagrangian or quasi-Lagrangian reference-frame where sweeping effects are 
absent || |4], |6|, Of course, from the experimental point of view, it is 
much harder to measure the velocity field in a Lagrangian frame than in the 
usual laboratory reference-frame. To our knowledge, the only results about 
multi-time velocity correlations are purely theoretical 0, |§ or numerical J|]. 
For this reason, in the following, we will limit to an experimental and theo- 
retical analysis of simultaneous -single time- multiscale correlation functions 
only. 

Multiscale correlations functions should play in turbulence the same role 
played by correlation functions in critical statistical phenomena. 
Recently, some theoretical work P, |6|, 0] and an exploratory experimental 
investigation || have been devoted to the behavior of single-time multiscale 
velocity correlations (MSVC): 

F p>q (r,R) = ((v (x + r) - v (x)) p (v (x + R) - v (x)) q ) 

= ((5 r v(x)r(5 R v(x)y) (3) 

with rd < r < R < L. When the smallest among the two scales r goes beyond 
the dissipative scales, r d , new properties of the correlation functions (||) may 
arise due to the non trivial physics of the dissipative cutoff [[[], |lOfl . From now 
on, we will mostly concentrate on correlation functions with both r and R 
in the inertial range, only in the last section we will address the important 
point concerned with the cross-over of the dissipative scale. Moreover, in 
order to simplify our discussion, we will confine our analysis to the case of 
longitudinal velocity differences. 

The main purpose of this paper is to review and to extend a recent experimen- 
tal and theoretical analysis of multiscale correlations (pf) |T5|. In particular, 



we present a systematic analysis of multiscale correlation functions in differ- 
ent experimental set-up, we also perform a critical comparison with the same 
observable measured in synthetic turbulent signals defined in terms of purely 
multiplicative random processes. 

The comparison with the synthetic signals will allow us to conclude that 
multiscale correlation functions are in quantitative agreement, with the pre- 
diction one obtains by using a pure uncorrelated multiplicative process for 
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the energy cascade, as long as both separations r, R are in the inertial range. 
As for the case when one of the two separations is already in the dissipative 
range we will critically review the two most important different predictions 
one can obtain imposing the dissipative cut-off using either multifractals 



or the GESS phenomenology HlQ| . Unfortunately, the actual state-of-the- 
art experimental dissipative-scales data does not allow to clearly distinguish 



among the two predictions [[Li], [12]. 



The paper is organized as follows. In section 2 we briefly review the ansatz 
that one simply obtains for MSVC (f|) by using a multiplicative random pro- 
cess for the inertial-range energy cascade. In section 3 we discuss sub-leading 
corrections induced by geometrical constraints which necessarily affects any 
MSVC for finite separation of scales. These geometrical constraints intro- 
duce sub-leading power laws behavior which may strongly interfere with the 
leading multiplicative predictions for finite separation of scales r/R ~ 0(1). 
In section 3 we also present our experimental data-analysis and the compar- 
ison with the synthetic multifractal field. In section 4 we briefly address the 
problem of dissipative correlation functions. Conclusions follow in section 5. 



2 Multiplicative random processes 

Stochastic cascade processes are simple and well known useful tools to de- 
scribe the leading phenomenology of the intermittent energy transfer in the 
inertial range. Both anomalous scaling exponents and viscous effects [[[J [10] 
can be reproduced by choosing a suitable random process for the multiplier, 
W(r, R), which connects velocity fluctuations at two different scales, R> r. 

The main idea turns around the hypothesis that small scale statistics 
is fully determined by a cascade process conditioned to some large scale 
configuration: 

5 r v(x) = W(r,R) ■ 5 R v(x) (4) 

where, requiring homogeneity along the cascade process, the random func- 
tion W should depend only on the ratio r/R. Structure functions are then 
described in terms of the W process: 

S p {r)=C p ([W{r/L)f) 1 
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where C p = ((5lv(x)) p ) if the stochastic multiplier may be considered almost 
uncorrelated with the large-scale velocity field. Pure power laws arise in the 
high Reynolds regime: in this limit we must have 



W 



r \ C(p) 

r) 



In the same framework, it is straightforward to give the leading prediction 
for the multiscale correlation functions (131): 



F P>q (r,R) 
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(5) 



which becomes in the hypothesis of negligible correlations among multipliers: 



F P , q {r,R) 




p+q\ 



S P (R) 



Sp+q(R) 



(6) 



This expression was for the first time proposed in || and later examined in 
more details in || where it was named "fusion rules". In the same article 
the authors proved that the fusion rule prediction gives the leading behavior 
of (0) when r / R — > as long as some weak hypothesis of scaling invariance 
and of universality of scaling exponents in Navier-Stokes equations hold. The 
name fusion-rules refers, probably, to the fact that thanks to the -supposed- 
uncorrelated nature of the cascade process every multi-scale correlation can 
be written in terms of single scale correlations -i.e. structure functions-. 
Let us notice that, beside any rigorous claim, expression (|j) is also the zero-th 
order prediction starting from any multiplicative uncorrelated random cas- 
cade satisfying S P (R) = (jW(j^)} p ^ S p (r). Let us also stress that the fusion 
rules prediction as stated in @ does not necessarily requires any scaling 
property of the underlying structure functions, a fact which suggests that 
the validity of the statement should be almost Reynolds independent. 



In this paper we want to address three main questions: 

• (i) whether the prediction @ gives the correct leading behavior in the 
limit of large scales separation: r/R ~ 0, 
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• (ii) if this is the case, what one can say about sub-leading behavior for 
separation r/R~ 0(1), 

• (iii) what happens to those observable for which the "multiplicative 
prediction" (|6|) is incorrect because of symmetry reasons. 

The last item comes from the observation that for correlation like : 

F 1 , q (r,R) = ((5 r v)(5 R v) q ) (7) 
the multiplicative prediction gives: 

Fi, q (r,R) = ^-S 1+q (R). 

Such a prediction is wrong because, if homogeneity can be assumed, Si(r) = 
for all scales r. In this case prediction @ does not represent the leading 
contribution. 

In the following we propose a systematic investigation of @ in high 



Reynolds number experiments [16, 17 , numerical simulation |18[ and syn- 



thetic signals |13| . The main purpose consists in probing whether multiscale 
correlation functions may show new dynamical properties (if any) which are 
not taken into account by the standard simple multiplicative models for the 
energy transfer. 



3 Data analysis 

In this section we present our main contributions by discussing the three 
items listed in the previous section and by presenting a detailed data-analysis 
in experiments at different Reynolds numbers, in numerical simulations and 
in synthetic signals. 

Experimental data sets come mainly from two different laboratories. We 
have analyzed data obtained in a wind tunnel (Modane) with Re\ = 2000, 
the integral scale was L ~ 20 m and the dissipative scale was = 0.31 mm. 
The second data set comes from a recirculating wind tunnel (ENS de Lyon) 
with a working section 3 m long and 50x50 cm cross section. Re\ involved 
in experiments were 400 (wake behind a cilinder) and 800 (jet turbulence). 
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Integral scales were 0.1m and 0.125m respectively whereas the dissipative 
scales were 0.15mm and 0.1mm. 

Synthetic signals are built in terms of a Wavelet decomposition with co- 



efficients defined by a pure uncorrelated random multiplicative process [[13 
In the following, the comparison between the synthetic field and the experi- 
mental data will play a central role in our discussion. Therefore, in Appendix 
A we briefly recall how a multiaffine field may be synthesized -and analyzed- 
in terms of a wavelet representation. 

In appendix A we prove that such a signal shows the fusion rules predic- 
tion (J^) and therefore it will turn out to be an useful tool for testing how 
much deviations from observed in experiments or numerical simulations, 
are due to important dynamical effects or only to unavoidable geometrical 
corrections. 

Let us proceed with a simple but basic observation. 
For any notice that for any 1-dimensional string of number (such as the 
typical outcome of laboratory experiments in turbulence) the multiscale cor- 
relations (|3]) feel unavoidable strong geometrical constraints. In particular, 
for any MSVC, with two velocities at the same spatial point v (x) and the 
two other velocities in a collinear geometry at spatial locations v(x + r) and 
v(x + R), like those analyzed in the following, we will always write down 
what we like to call the "Ward-Identities" (WI): 

S p (R-r) = ([(v(x + R) - v(x)) - (v(x + r) - v(x))] p ) (8) 
= £ b(k,p)(-) k F k , p _ k (r,R), (9) 

k=0,p 

where b(k,p) = p\/[k\(p — k)\]. 
For example, for p = 2 we have 

2F ltl (r,R) = S 2 (r)+S 2 (R)-S 2 (R-r)~ 
r \ £( 2 ) / r \ 

_) +0(-)\.S 2 (R) (10, 

where the latter expression has been obtained by expanding S 2 (R — r) in the 
limit r/R — > 0. 
For p = 3 we have 

S 3 (R -r) = S 3 (R) - S 3 (r) + 3F 2>1 (r, R) - 3F 1>2 (r, R) 
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The "Ward-Identities" will turn out to be useful for understanding sub- 
leading predictions to the multiplicative cascade process. One may argue that 
in geometrical set-up different from the one specified in @ the same kind 
of constraint will appear with eventually different weights among different 
terms. 

The most important result one must extract from (||) is that the mul- 
tiscale correlation functions, as stated in (JJ), may not be a perfect scaling 
functions even in the limit of very high Reynolds number. Indeed, the WI 
tell us that MSVC with different order of velocity moments must be con- 
nected unavoidably one with the other: which would be in contrast with the 
assumption that each MSVC should be determined by a single power law 
behavior. 

The main result presented in this work is that all multiscale correlations 
functions are well reproduced in their leading term, — > 0, by a simple 
uncorrelated random cascade ([]) and that their sub-leading contribution, 
~ 0(1), are fully captured by the geometrical constrained previously dis- 
cussed, namely the "Ward-Identities". 

The recipe for calculating multiscale correlations is the following: first, 
apply the multiplicative guess for the leading contribution and look for geo- 
metrical constraints in order to find out sub-leading terms. 
Second, in all cases where the leading multiplicative contribution vanishes be- 
cause of underlying symmetries, look directly for the geometrical constraints 
and find out what is the leading contribution applying the multiplicative 
random approximation to all, non-vanishing, terms in the WI. 



3.1 Fusion rules: even moments 

Let us check the fusion rules prediction (|6|) for even moments p, q — 2, 4, 

In order to better highlight the scaling properties we will often use in the 
following, F Ptq (r, R), the MSVC compensated with the fusion- rule prediction: 

P , r m_ F r ,,(r, R) ■ S„(R) 



In order to compare experiments with different Reynolds numbers we may 

R-r 
L-r ' 



use as independent variable in our plot the quantity: x(R) = j- 1 :, where with 
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L we intend the integral scale of each different experiment. In this way, by 
fixing the small scale r = 5r] and by changing r < R < L for each set of 
data we have a variation of < x(R) < 1. In Figure (1) we have checked 
the large scale dependency by plotting the compensated MSVC as a function 
of x = at fixed, r, for p = 2, q = 2 and different Reynolds numbers 
(experimental and numerical). 

The expression (||) predicts the existence of a plateau (independent of R) at 
all scales r < R < L where the leading multiplicative description is correct. 

From Figure (1) one can see that experiments with low Re numbers show 
a slightly poor plateau. In particular the DNS (Direct Numerical Simulation) 
with Re\ ~ 40 does not show any plateau. The absence of a plateau is con- 
nected to the overwhelming geometrical effects present at such low Reynolds 
numbers (see below). 

For this reason, in the following figures we will show only experimental data 
from Modane wind tunnel, which have the highest Re number we can access. 

In Figure (2) one can see the compensated correlation functions for two 
different set of moments. In the limit of large separation R — > L at fixed r, 
we indeed see a tendency toward a plateau. On the other hand, there are 
clear deviations for r/R ~ 0(1). 

The same behavior is seen in Figure (3) for the same compensated quantities 
fixing the large scale R and by changing the small scale r. 
Such deviations show a very slow decay as a function of the scale separa- 
tion. The decaying is so slow that a clear plateau is seen only for the largest 
Reynolds number available. The question whether the observed finite-size 
corrections have an important physical origin or not is therefore of primary 
importance. 

In order to understand the physical meaning of the observed deviations 
to the fusion rules (0), we compare, in Figure (4), the experimental data 
against the equivalent quantities measured by using the synthetic signal. 
We notice an almost perfect superposition of the two data sets, indicating 
that the deviations observed in real data can hardly be considered a "dy- 
namical effect" . 

Using the WI plus the multiplicative ansatz for the leading behavior of all 
correlation functions for p = 4 we quickly read that the leading contribution 
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to F 2j 2 is 0(r^ 2 )) • 0(i?^ 4 ) ^ 2 ^), while sub-leading terms scale as 0(r c ^), 
and as 0(r« 3 )) • 0(R^~^). 

This superposition of power laws is responsible for the slowly-decaying cor- 
relations in Figures (1-4). 

In Table 1 we summarize the leading and sub-leading contributions that may 
be inferred from the WI for the standard MSVC with p = 2, q = 2. Similar 
tables can be constructed for any other even MSVC. 

The result so far obtained, i.e. that both the experimental data and the 
synthetic signal show the same quantitative behavior, is a strong indication 
that multiscale correlation functions, at least for even order of the moments 
-i.e. in all cases where the signal is not affected by cancellation problems-, 
are in good agreement with the random multiplicative model for the energy 
transfer. 

An even stronger prove of this statement comes from the analysis of multi- 
scale correlations in terms of the coefficients obtained by a Wavelet analysis 
of the experimental signal (see appendix A). 

The Wavelet coefficient a,^ may be seen as the representative of a velocity 
fluctuation at scale r = 2~ j and centered in one of the k = 1, 2, • • • , 2 j spatial 
point chosen equispaced in the original total length of the signal. 
With this interpretation in mind, we may think at the Wavelet coefficients as 
the ideal observable which minimize the geometrical constraints and therefore 
as the ideal cases where one can test the idea that behind the multiplicative 
process there are only geometrical constraints. In other words, in terms of 
the coefficients obtained by a Wavelet analysis of the experimental signal, the 
multiscale correlation function should show the fusion rules prediction for a 
range of scales much wider than for the velocity increments, i.e. geometri- 
cal constraints, which introduce sub-leading power-laws decaying, should be 
minimized. 

In Figure (5) we indeed prove that this is exactly what happens. In Fig- 
ure (5) we show the equivalent of F 2;2 (r, i?) built in terms of the Wavelet 
coefficients: 

F 2 y(r = 2~*,R = 2"'") = (\a hk \ 2 \a r , k ,\ 2 ) (12) 

In the figure we plot, as in the previous figures, the compensated correlation, 
obtained from the Wavelet coefficients, at fixed small scale and at changing 
the large scale. In Figure (6) the same quantities are plotted at changing the 
small scale. As it is evident, the finite scale-separation effects visible in the 
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standard MSVC have here disappeared. 

The plateau is reached immediately after, say, one fragmentation step. 



3.2 A case where fusion rules fail 



For multiscale correlations where the direct application of the random-cascade 
prediction is useless -because of the translation symmetry-, like Fi >g (r, R), we 
suppose that the main leading contribution is simply due to the geometrical 
constraints. In other words, we say that as soon as the main leading effect in- 
duced by the presence of a multiplicative random energy transfer is depleted 
because of symmetry reasons, the sub-leading contributions induced by the 
geometry becomes the leading contributions. 

In order to give a prediction for such class of MSVC we therefore use the WI 
applying the multiplicative prediction to all terms, except the F\ >q . 
One obtains the expansion: 



Fi >q (r,R) 



r \ / r \ 



r \ C(4) 



r \ C(9+i)' 



(13) 



which coincides when q = 1 with the exact result (10) using ((3) = 1. 

In Figure (7) we show the experimentally measured Fi^ and the fit that 
we obtain by keeping only the first two terms of the expansion in fll3|). The fit 
has been performed by imposing the value for the scaling exponents C(2), C(3) 
measured on the structure functions, i. e. only the coefficients in front of the 
power laws have been fitted. As one can notice, the fit works perfectly in the 
inertial range. Let us remark that the correlation changes sign in the middle 
of the inertial range, which is a clear indication that a single power-law fit 
(neglecting sub-leading terms) would completely miss the correct behavior. 

Next we consider the WI for p = 3. Due to the fact that Ss(r) ~ r in the 
inertial range, one can easily show that the WI enforces F\% ~ F21. Therefore 
we can safely state that also correlation functions of the form F Pj i feel non 
trivial dependency from the large scale R, proving that the prediction given 
in using isotropy arguments is wrong. 
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3.3 Fusion rules: odd moments 



For the most general MSVC involving odd moments of velocity increments, 
F Pj g(r, R) with p, q — 3, 5, 7, . . .. the situation is slightly more confused. The 
problem stems from the fact that the fusion rules contribution to this correla- 
tion feels indeed the skewed part of the process which is order of magnitudes 
less important than the even part. For example, the multiplicative contri- 
bution to the correlation Fs t i(r,R) would be Ss(r) / Sz(R) ■ S^R) which is 
different form zero only due to the fact that the process for the longitudinal 
velocity correlation is skewed. 

The weakness of the signal from the multiplicative contributions makes these 
class of correlation functions very hard to analyze from the point of view of 
scaling. Here, the geometrical constraints may well be more important, in a 
large range of scale separation, than the fusion rules prediction. For example 
in Figure (8) we plot the standard MSVC for p = 3, q = 1 and the same 
correlations but with moduli of velocity increments, such as to get rid, in 
the second case, of cancellation effects. As it is evident, the two correlation 
have a very different amplitude as soon as the scale separation becomes im- 
portant and it is hard to say whether the MSVC without moduli follow the 
fusion rules prediction for large scale separation or not. On the other hand, 
the correlation with absolute values does follow the multiplicative prediction 
reaching a plateau after the usual finite size transient as the ordinary even- 



A high statistics and high Reynolds number investigation of such a class 
of correlation may well be of some interest in order to elucidate whether the 
odd part of velocity increments follows the same physics of the even part or 
not. 

4 Dissipative Physics 

In this section we discuss the application of fusion rules in the dissipative 
range. We will be mainly interested in the following two quantities: 



MSVC. 



MR) 



(Av(x) ■ Sv R {x) n ) 



(14) 




(15) 
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where Av(x) is the laplacian computed at the point x, Svr(x) = v(x + R) — 
v(x) and T(x) is the velocity gradient computed at x. In order to simplify 
the discussion we restrict to the one dimensional case, namely the Laplacian 
and the gradient are computed in one dimension and velocity differences are 
longitudinal. Moreover we restrict our analysis to the cases of n odd and 
n +p even. Our findings will anyway be valid in the most general case. The 



scaling properties of A n and B p n have been investigated in ||, |TT[ and |[L5f| . 

We start by considering the scaling properties of A n . By its definition we 
have: 



a ,n\ i- / f v(x + r) +v(x — r) —2v(x)\ , , . .. n \ 

A n {R) = Jim((- y — 2 ^1 • (v(x + R) -v(x)) n ) = 

= \imr' 2 (F ljn (r,R) + F hn (-r,R)) (16) 

In order to understand how equation (|T^) works, we compute the easiest 
quantities, i.e. A\ and A 3 . By using equation (|10D, we obtain: 

_2 

Ai(R) = hm T —(2S 2 (R) + 2S 2 (r)-S 2 (R-r)-S 2 (R + r))= (17) 

r— >0 2 

,-2 (a /_\ l 2 d 2 S 2 (R) , ?> 



limr- 2 [S 2 (r) - if-^ + (r*) j = (18) 



In equation (|l9|) , we have used the relation: 

The computation of A 3 is similar and we find, using (§): 

A 3 (i?) = limr- 2 ^(25 4 ( J R) + 25 4 (r) + 12F 2i2 (r, i2) - 4F 3il (r JJ R) + 
-4F 3 ,i(-r, i?) - S 4 {R - r) - 5 4 (i2 + r)) = 

= limr- a (3F 2 ^(r,J2)- V^^ + ofr 3 
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In equation (|2T| ) we used the definition of -82,2, namely 



£2,2 (R) = lim r- 2 F 2 , 2 (r, R) (22) 

r — ^0 

At this point it is quite easy to find the most general expression for A n , which 
is 

MR)=nB^ W--L-iS$* (23) 

Equation fl23|) is an exact results which is independent on any physical as- 
sumption on the fusion rules. The last term on the r.h.s. of (|23|) becomes 
small for R in the inertial range. On the other hand, for small value of R, 
i.e. for R — > 0, the last term of the r.h.s. of (|23|) cannot be neglected. In 
particular, for R — ■> 0, an explicit computation, either using (|23|) or flll]) , 
gives -after cancellations of leading terms in the RHS of (|23|)-: 

A n {R) ~ O (R n+1 ) (24) 

In order to complete our computation for A n , we need an estimate for B 2n . 
There are in principle two ways to compute B vn : the first one using the 
multiscaling approach |L4]], the second one using the GESS theory discussed 

m. 



in 



We first analyze the case of multiscaling. In this case, one can use the 



approach of multiplicative processes with multiscaling viscous cutoff [[14 
Namely, for the correlation B 2in {R) = \ {.d x v) 2 (5rv(x))) one obtains: 



B2AR) ~ ^Rv(x)) n ^ (25) 

where r d is the dissipative scale. In the multifractal interpretation we assume: 
5 Td v = (r d /R) h ■ 5rv with probability Ph(r d ,R) = (r d / R) 3 ~ D ^ h \ Following 
[ |i~4H we have: 



S u v ■ r d ~ ( — j 5 R v -r d ~v. (26) 
Inserting the last expression in the definition of B 2 ^ n (R), we finally have: 

fl2ww (^.(_i_) « „w* (27) 
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where we have used the fact that the multifractal process is such that 
v ((d x v) 2 ) — > 0(1) in the limit v — > 0. Expression ( p7|) coincides with the 
prediction given in [0] . The above computation are easily generalized for any 
{(d x v) p (5 R v(x)) g ). By using (p7j) and (p3|) we finally obtain: 

MR) = nC n ^§-^-^^ (28) 

Let us note that, for R — > 0, equation (p8j) predicts that A n (R) ~ O (-R n_1 ) 
which violates equation ([Zip. 

We now compute A n (i?) by using the GESS approach discussed in [[K] 



In this case the computation of B 2)n can be easily done by noting that, within 
the GESS approach, the fluctuations of the dissipation scale are confined in 
the range where 5rv ~ R. This implies that, for what concerns the scaling 
properties of £?2,n(-R), the effect of a fluctuating dissipation scale can be 



disregarded. Following []10| , after a long but straightforward computation, 
we obtain: 

Sa " {R) = Dn {T2) sXR) R) (29) 
where T is the velocity gradient and D n is a constant. Equation (^) can be 
easily understood by noting that, within the GESS approach, -F 2in (r, R) ~ 

^^sfcR) 2 ^ ^ or an y v& l ues °f r ari d ^> *.e. a ^ so i n ^e limit r — > 0. Using (|^) 
we finally obtain: 

MR) = nDn S 2 (R) - n + 1 d& m 

For R -> 0, using the estimate 5 n (i?) ~ (T n ) i? n + (i? n+2 ), and the fact 
that D n = 1 for R = 0, we can easily show that equation (PD| ) satisfies the 
constrain (p4j). 

From an experimental point of view, it is extremely difficult to distinguish 
between the two predictions ( p8|) or (|30|). We note that the experimental 
and numerical analysis discussed in [[□], |I2|| , has been done neglecting the 
second term on the r.h.s. of (^). Also, the experimental analysis performed 
T0[1 seems to indicate that multiscaling effects are not observed in real 



in 



turbulence. At any rate, no definitive conclusions can be drawn from existing 
experimental data. 
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5 Conclusions 



Let us summarize what is the framework we have found until now. 
Whenever the simple scaling ansatz based on the uncorrelated multiplicative 
process is not prevented by symmetry arguments, the multi-scale correla- 
tions are in good asymptotic agreement with the fusion rules prediction even 
if strong corrections due to sub-leading terms are seen for small-scale sep- 
aration r/R ~ 0(1). Sub leading terms are strongly connected to the WI 
previously discussed, i.e. to geometrical constraints. In the other cases (i.e. 
Fi tq (r,R)) the geometry fully determines both leading and sub-leading scal- 
ing. 

All these findings, led us to the conclusions that multiscale correlations func- 
tions measured in turbulence are fully consistent with a multiplicative, almost 
uncorrelated, random process. 

Nevertheless, the strong and slowly-decaying sub-leading corrections to the 
naive multiplicative fusion rules predictions are particularly annoying for any 
attempts to attack analytically the equation of motion for structure func- 
tions; in that case, multiscale correlations at almost coinciding scales are 
certainly the dominant contributions in the non-linear part of the equations 
M. Indeed, as shown in an analytical calculation for a dynamical toy model 
of random passive-scalar advection ||19|| , fusion rules are violated at small 
scale-separation and the violations are relevant for correctly evaluating the 
exact behavior of structure functions at all scales. 

Finally, let us remark that the standard multiplicative process may not 
be the end of the story, i. e. the dynamics may be more complex than what 
here summarized. 

For example, one cannot exclude that also sub-leading (with respect to the 
multiplicative ansatz) dynamical processes are acting in the energy transfer 
from large to small scales. These dynamical corrections must be either neg- 
ligible with respect to the geometrical constraints or, at best, of the same 
order. Also, as shown in this paper, the odd correlation functions are not jet 
under control: higher Reynolds number experiments, with higher statistics, 
are needed. 

The question connected to the transfer properties of quantities "orthog- 



onal" to the energy may reveal different physical mechanisms [23|. What 



happens for all those multiscale correlation functions which feel a non-trivial 
helicity dependency for non-parity invariant flows is in this framework an 
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open question. 

For what concerns fusion rules involving velocity gradients or laplacian 
and velocity differences, we observe that there are controversial arguments 
leading to different predictions. It is difficult to distinguish which predic- 
tions is really observed in real turbulence, because experimental data at 
large Reynolds number do no resolve the far dissipative range with enough 
accuracy. 

Other possible candidates to investigate the previous problems are shell 
models for turbulence, where geometrical constraints do not affect the energy 
cascade mechanism. 

We acknowledge useful discussions with A.L. Fairhall, V. L'vov and I. 
Procaccia. M. Pasqui is kindly acknowledged for his help in the analysis of 
the synthetic signal. We are indebted to Y. Gagne for having allowed us 
the access to the experimental data. L.B. and F.T have been supported by 
INFM (PRA TURBO). G. R. C. and S. C. acknowledge support by ECOS 
comitee and CONACYT under project M96-E03. 



6 Appendix A 

We build up a 1-dimensional synthetic signal according to a random mul- 
tiplicative process defined in a dyadic hierarchical structure as originally 
introduced in |0| (for a review and references see also [pOfl ). 
Let us consider a wavelet decomposition of the function 4>(x): 

oo 

0( x ) = a jik i/; j>k (x) (31) 

j,k=0 

where ipj^{x) = 2- J ' 2 '0(2 J x — k) and ^(x) is any wavelet with zero mean. The 
above decomposition defines the signal as a dyadic superposition of basic 
fluctuations with different characteristic widths (controlled by the index j) 
and centered in different spatial points (controlled by the index k). For 
functions defined on N = 2 n points in the interval [0, 1] the sums in (|31|) are 
restricted from zero to n — 1 for the index j and from zero to 2 J — 1 for k. 



In it has been shown that the statistical behavior of signal increments: 

< \5<f){r)\ p >=< \<j)(x + r) - <f){x)\ p >~ 
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is controlled by the coefficients ct^. By defining the a coefficients in terms 
of a multiplicative random process on the dyadic tree it is possible to give an 
explicit expression for the scaling exponents C(p)- F° r example, it is possible 
to recover the standard anomalous scaling by defining the a's tree in term of 
the realizations of a random variable rj with a probability distribution P(r)): 



"0,0 



"i,o = ??i,o«o,o; = ii,i«o,o; 

"2,0 = ?72,0 «1,0; "2,1 = V2,l «1,0; "2,2 = %2«l,i; «2,3 = *,3«1,1, (32) 

and so on. Let us note that in the previous multiplicative process different 
scales are characterized by different values of the index j, i.e. rj = 2~ J . If 
the rjj^ are indipendent identically distributed random variable it is straight- 
forward to realize that are random variables with moments given by: 

< \ a . k \v >= r-^tf) =7 .C(p) (33) 

where the "mother eddy' ao,o has been chosen equal to one. In (^) with 7TT 
we intend averaging over the P(rj) distribution. In [JT3] it has been shown 
that also the signal <j)(x) has the same anomalous scaling of ([33]). 
The same arguments used in order to prove that the field <j)(x) has an anoma- 
lous scaling can be invoked to show also that the fusion-rules prediction flfj) 
are satisfied -at least for large scale separation-. 

On the other hand, it is a trivial matter to realize that the above signal will 
show exactly, and for any separation of scale, the fusion-rules prediction if 
expressed for the wavelet coefficients atj^. For example, let us consider two 
wavelet coefficients at different scales rj < ry and let us chose the k-s indices 
such that the two coefficients refer to two spatially overlapping wavelets, then 
it is trivial to realize that, due to the multiplicative nature of the wavelet 
coefficients, we have: 




<h>f l«^H= -M 4™ (34) 



which shows that the fusion rules prediction is satisfied exactly for any sep- 
aration of scales as long as the two fluctuations are chosen with overlapping 
distances. In the case the two distances are not overlapping, deviations from 



18 



the fusion rules prediction are certainly seen in the synthetic field due to 
the dyadic -ultrametric- nature of the underlying structure. The question 
whether such deviations may be seen also in the experimental data is an 
interesting point which is outside the scope of this paper (see for similar 
problems 
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FIGURES and TABLE CAPTIONS 

TABLE 1: 

Leading (first column) and sub-leading (second column) contribution to the 
different multi-scale velocity correlations entering in the WI written for p — 4. 
Notice that all the leading behaviors have been obtained by using the multi- 
plicative ansatz (when applicable). The sub-leading behaviors are consistent 
with the constraints imposed by the WI. 

FIGURE 1: 

Compensated MSVC F 2 ^(r, R) at fixed r and changing x(R) = for differ- 
ent experiments and numerical simulation: ( x ) Direct Numerical Simulation 
(Re x = 40), (+) Jet (Re x = 800), (*) Modane (Re x = 2000), (□) Wake 
(Re x = 400). 

FIGURE 2: 

Compensated MSVC F PjQ (r, R) at fixed r and changing the large scale R for 
p = 2,q = 2 (+) and p = 2, q = 4 (x). 

FIGURE 3: 

Experimental compensated MSVC F p q (r, R) / S p+q (R) ■ S p (r) at fixed R and 
changing the small scale 1/r for p = 2, q = 2 (+) and p = 2, q = 4 (x). 

FIGURE 4: 

Comparison between experimental and synthetic compensated MSVC, F Pjq (r, R) 
at fixed r and changing the large scale R for p = 2, q = 2: (+) synthetic and 
(x) experimental. For p = 2, q = 4: (*) synthetic and (□) experimental. 

FIGURE 5: 

Comparison between real space (+) and wavelet analysis (x) of the exper- 
imental data set from Modane. Compensated F 2 ,2 is shown for fixed r at 
varying R. 

FIGURE 6: 

Comparison between real space (+) and wavelet analysis (x) of the exper- 
imental data set from Modane. Compensated F 2i2 is shown for fixed R at 
varying r. 
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FIGURE 7: 

Experimental Fi t 2(r,R) at fixed r = 16 and at varying R. The integral 
scale L ~ 1 x 10 4 r^. Let us remark that the observed change of sign in the 
correlation implies the presence of at least two power laws. The continuous 
line is the fit in the region r < R < L obtained by using only the first two 
terms in ([13]). 

FIGURE 8: 

Comparison between compensated F 3j i odd MSVC with absolute values (+) 
and without ( x ) . Data are shown for fixed r at varying R. It is evident that 
the odd MSVC with absolute values has the same behavior of even MSVC, 
while the one without absolute value does not follow the same behavior. 
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Figure 2: Benzi et al. 
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Figure 6: Benzi et al. 




Figure 7: Benzi et al. 
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